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(w) A method and system for ascertaining density features of underground earth regions, and a detector module usable 
therein for detecting the flux and trajectories of muons. 

© A plurality of detectors (12) is disposed within s bore 
hole for escertfiintng the frequency and trajectory of high- 
energy muons that penetrate through the earth. This infor- 
mation Is converted by data processor (42) into an integrated 
density value of the earth for incremental angular lines of 
sight extending through spaced-apart elevations along the 
detectors (12). The data processor (42) uses a tomographic 
analysis process to combine together alt of those integrated 
density values (or all of the elevations of interest to calculate 
the densities of discrete spatial volumes composing the 
underground region being investigated. Detectors (12) pro- 
duce electrical signals indicative of the trajectory of the 
passing muon. The signals are recorded by tirne-to-dighal 
oonverter l ,r TDC") units (30). if the signals appear to have 
possibly originated from a high-energy muon, a micro- 
processor (381 retrieves the information collected by TDC 
units (30) end transmits it toe data logger (40) located at the 
earth's surface. The data processor (42) then utilizes tomo- 
graphic analysis to make the localized density determina- 
tions discussed above. 
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A METHOD AND SYSTEM FOR ASCERTAINING DENSITY 
FEATURES OF UNDERGROUND EARTH REGIONS, AND A 
DETECTOR MODULE USABLE THEREIN FOR DETECTING 
THE FLUX AND TRAJECTORIES OF MUONS 

Field of the Invention 
The present invention relates to methods and apparatus for 
investigating the characteristics of underground earth formations by determining 
underground density patterns, and more particularly to using cosmic ray muons to 
ascertain mean or integrated densities along numerous trajectories, and 
subsequently using tomographic analysis to convert these mean densities into 
local density values over the entire region covered by the muon trajectories. 

Background Art 

In the past, characteristics of underground earth formations have 
been investigated by use of various techniques that can be divided into those 
which require drilling and those which do not. In the latter category are such 
surface methods as gravimetric surveys that detect the earth's gravitational 
anomalies and seismic surveys that utilize shock waves to investigate changes in 
the eartlVs density. While these methods are useful, they tend to exhibit large 
scale features and do not yield detailed localized information. In the former 
category, drilling techniques through, for instance, coring or mud logging, yield 
abundant information, but are highly, localized covering a small cylindrical region 
surrounding the hole. Thus, a combination of techniques may be required; for 
example, the existence of an ore body may be determined by gravimetric 
measurements, but its precise delineation may require the drilling of numerous 
holes. 

Several other known down-hole study methods use gamma rays, 
whether naturally occurring or induced by neutron activation. Because gamma 
rays have ranges of a few meters in earth, the methods depending on them yield 
information only on a region of the same size. Examples of gamma ray detection 
systems designed to determine the composition of underground earth formations 
and distinguish between rock, sand, silt, water and other components include U.S. 
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Patents 2,332,873; 2,984,745; and 2,998,521. Other systems are specifically 
intended to provide information regarding the porosity of the strata through 
which a bore hole or well extends. Examples of such systems are found in U.S. 
Patents 3,509,342; 4,005,290; and 4,035,639, As disclosed in U.S. Patents 

5 2,963,587; 3,156,822; and 3,278,748, a further type of system utilizes a series of 
gamma radiation detectors to analyze the chemical elements that are present in 
the earth formations surrounding a well bore. 

Gamma ray detection systems of the nature discussed above 
employ various types of gamma radiation counting or detecting devices to 

10 measure both naturally occurring and induced gamma radiation. The use of 
Geiger-Muller detectors is disclosed by U.S. Patents 2,332,873; 2,963,587; 
2,998,521; 3,156,892; and 3,278,748. Scintillation counter detectors coupled with 
photomultipliers are employed in U.S. Patents 2,830,187; 2,963,587; 2,984,745; 
3,156,822; 3,278,748; and 3,509,342. Gamma ray detection systems incorporating 

15 ionization chamber detectors are set forth by U.S. Patents 2,332,873 and 
4,035,639. A gamma ray detection system including a proportional counter 
detector is disclosed by U.S. Patent 4,005,290. U.S. Patent 2,721,943 utilizes a 
combination scintillation/Geiger-Mulier detector composed of a Geiger-Muller 
type photoelectric detector housed within an inner cylinder, which is disposed 

20 within an outer cylinder containing liquid luminophor. The scintillations emitted 
by the liquid luminophor from Incident gamma radiation are amplified by the 
light detector and then converted into electrical pulses. 

Detectors utilized in gamma ray detection systems have been 
designed to monitor only gamma rays reaching the detector from particular 

25 directions. In U.S. Patents 2,830,187; 2,984,745; and 3,156,822, the detectors 
screen out all gamma radiation except radiation approaching from a direction 
generally perpendicular to the length of the hole. 

As stated above, gamma ray techniques are only capabl of 
providing information concerning conditions within a few meters of the hole, and 

30 more information can generally be obtained only by drilling more holes. In 
contrast, cosmic ray muons have a spectrum of energies that extends well beyond 
a million MeV, and such muons, especially the more energetic ones, have ranges 
in earth measured in kilometers. The study of the transmission of such muons 
can therefore yield information concerning the density of the earth over 

35 distances of this general magnitude* 

The sources, trajectories, isotropy, energy spectrum, composition 
and ther properties of cosmic ray muons have been exhaustively studied in an 
academic context. For instance, Crouch et al. in "Underground Measurements f 
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Cosmic Ray Muons and Neutrinos", Proceedings of the Llth International 
Conference on Cosmic Rays , Budapest 1969, Paper MU-21, r ports installing 
scintillation detectors In deep underground locations to investigate the trajec- 
tories at which muons and neutrinos strike the earth, Keuf fel et al. in "Zenith- 
5 Angle Distribution of Ultra-High Energy Muons", Proceedings of the llth 
International Conference on Cosmic Rays, Budapest 1969, Paper MU-31, 
discusses placing muon detectors in underground locations to study the intensities 
of incident muons as a function of their zenith angle. 

As discussed by Alvarez in "A Method to 'X-Ray T the Egyptian 
10 Pyramids to Search for Unknown Chambers", Adventures in Experimental 
Physics , p. 175 (1972), cosmic ray muons have been utilized in attempt to search 
for hidden chambers in the Second Pyramid of Chephren in Egypt. The 
trajectories of the muons passing through the pyramid were monitored by two 
vertically spaced-apart spark chambers. Each spark chamber was composed of 
15 two flat metal plates disposed horizontally relative to each other and spaced 
approximately one centimeter apart. When a muon penetrated the two plates, a 
high voltage was applied between the plates causing a spark to jump from one 
plate to another along very nearly the same trajectory as the passing muon. A 
separate triggering system was required to apply the high voltage across the 
20 metal plates at the correct time. The locations of the sparks were recorded by 
two cameras aimed at the gap between the plates of the spark chambers with the 
cameras disposed at right angles to each other. In a later version of the 
apparatus, the coordinates of the spark were recorded directly on magnetic tape 
by a nickel wire magnetostrictive readout. 
25 The present invention involves determining the mean or integrated 

density of the earth along paths that muons travel through the earth and then 
using tomographic analysis to convert these mean densities into local density 
values. In the past, tomographic techniques have been used in medical 
applications. In the medical context this technique employs the principle that X- 
30 rays, traversing the body, are partially absorbed by the tissue in their path, and 
the transmission of these X-rays therefore indicates the total integrated density 
of tissue along the path of the X-rays. Ordinary X-ray procedure involves using a 
fan-shaped beam of X*-rays emanating from a fixed source. Such an exposure 
yields the integrated tissue density along all the lines traversed by the fan. 
35 Typically, in tomography, the sourc of the X-rays and the sensing plane, for 
instance a film or other detector, are rotated together while th object, a head 
or a whole body, remains fixed. As a result, the integrated density along lines 
having many different directions is obtain d. A mathematical procedure is then 
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employed to convert the totality of such integrated measurements into localized 
density values* Once this conversion has been performed, the local densities thus 
obtained can be displayed in many different forms, one of which is a contour map 
showing local densities in any chosen plane. 
5 In the present invention, cosmic ray muons are used in place of X- " 

rays since their range is appropriate to the size of the subjects being studied. 
Cosmic ray muons strike the earth from all directions, and as long as the * 
direction from which each muon arrives is properly measured, this isotropy plays 
the same part as the rotation of the source of X-rays in medical tomography* 
10 The data thus obtained yields the integrated densities along lines arriving at a 
particular detector from a range of directions. An array of detectors may be 
arranged, for instance, to extend down a hole instead of in a plane as in a 
medical situation. But by a substantially similar mathematical procedure, the 
totality of line densities thus obtained from all of the detectors can be converted 
15 into detailed local densities. This information can then be displayed in any 
convenient form, including the local densities in any plane. 

General discussions of tomography in a medical context appear in 
the articles "Image Reconstruction from Projections" by Richard Gordon, 
Scientific American , pages 56-68 (October 1975); and "Medical Imaging Sys- 
20 terns", Rolland Redington and Walter Berniger, Physics Today , pages 36~4<f 
(August 1981), Prior art patents generally pertaining to tomography systems and 
specifically concerning different tomographic analysis methods are disclosed by 
U*S. patents 4,044,240 (concerning a transverse section tomography system for 
providing a reconstructed image of a planar section of a subject); 4,072,289 
25 (pertaining to an analog method and apparatus for producing an axial tomo- 
graphic image of the human body); 4,135,247 (involving a tomography system for 
processing data from a plurality of projections in essentially real time to produce 
a reconstructed image immediately following completion of a scan); 4,149,249 
concerning a method and apparatus for tomographic analysis wherein the linear 
30 projections are obtained by passing radiation in quasi-planar form through the 
quasi-plane of the object slice at various angles and measuring the amounts of 
radiation not absorbed or scattered by the object slice); 4,168,435 (concerning a 
method for preprocessing radiological information by means of electronic analog 
techniques for subsequent use in a process for reconstructing a three-dimensional 
35 image of a desired laminal of a subject); 4,205,375 (involving a specific 
tomographic analysis procedure for reconstructing a three-dimensional image 
utilizing Fourier transformations); 4,212,062 (pertaining to a tomographic 
imaging process in which a convolution integration over an infinite region is 
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approximated by an integration over a finite region by employment of a 
predetermined weighing function); 4,217,6^1 (concerning an image-processing 
apparatus that makes corrections for polychromatic distortion and images 
produced by the nonlinear interaction of body constituents with a polychromatic 
X-ray beam); 4,222,104 (relating to a process for improving the accuracy of the 
reconstruction portion of a tomographic process by first processing the X-ray 
attenuation data using a procedure and electronic circuits for fan-distributed 
radiation paths and then reprocessing the data using the same algorithm 
procedure and circuits for parallel radiation paths); 4,233,622 (concerning a 
10 corrected layergram technique); 4,259,725 (pertaining to a system for generating 
cursors that overlie and move on a preliminary computerized tomography 
generated X-ray view that is already on a display screen); and, 4,293,912 
(concerning an apparatus and method for performing a convolution operation on 
at least a part of the sets of electronic data which represent absorption in the 
15 fan-shaped regions without binning the electronic data into parallel ray data sets 
with a convolution function to form convoluted data). 

Disclosure of the Invention 
The present * invention relates to a method and an apparatus for 
utilizing cosmic ray muons to ascertain the local density features of large 
20 underground earth formations to assist, for instance, in locating mineral deposits, 
promising oil or natural gas reservoirs, coal horizons, voids, structural anomalies 
and other features that can be related to local density variations. Cosmic ray 
muons continuously strike and penetrate the earth from all directions. These 
muons possess the falling energy spectrum. Like other charged subatomic 
25 particles, muons lose energy by ionization processes. The energies of muons 
diminish as they travel through the earth, and more specifically, the rate at 
which muons lose their energy, arid thus their attenuation, is related to the type 
and density of the material through which they pass. For a given energy, muons 
can travel further through less dense material than through relatively more dense 
30 material. Thus, it is expected that less attenuation of the muon flux will occur 
in portions of an underground earth formation that is less dense than in sections 
which are of higher density, since lower energy muons are also capable of 
penetrating the less dense portions of the earth formation. In other words, for a 
given time interval, more muons will arrive at a detector alon" a path of lower 
35 density than along a path of higher density. Thus, by monitoring the frequency 
and traj ctories at which mu ns pass through the arth, the density features of 
the underground earth formation may be explored. These concepts are used to 
investigate the local density characteristics of large volume underground arth 
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formations, 

A plurality of muon detectors are placed at underground locations 
to monitor the flux and trajectories of incident muons striking the detectors. 
From this data the integrated density along incremental angularly separated lines 
5 of sight extending through spaced-apart elevations or elevational ranges of the 
detector can be determined. For each individual elevation or elevational range 
along the detectors, the integrated density along a particular direction Is the 
result of the densities of all of the different local density variations along the 
path. 

10 By the present invention the size and location of individual 

structural anomalies or other density-related features may be determined. This 
is accomplished by using a mathematical procedure to convert the totality of all 
of the integrated density measurements along all of the lines of sight for all of 
the selected elevations along the detectors into localized values. In other words, 

15 the density at each point or spatial volume comprising the region under 
investigation is determined. These localized density values can be displayed in 
various manners, such as a contour map showing local densities upon any chosen 
plane in a manner similar to that commonly generated in medical examinations. 

According to the method of the present invention, each time a 

20 muon penetrates the earth and passes through a detector, a first set of electrical 
signals is produced on particular sensing wires disposed within the detector and a 
second set of electrical signals is induced on certain pickup wires also located 
within the detector. The relative timing of the individual signals of the first set 
and the particular wires in which the signals are produced are indicative of the 

25 azimuth angle of the muon; and the timing of the individual signals of the second 
set and the particular pickup wires in which the signals are induced are indicative 
of the zenith angle of the muon and the elevation along the height of the 
detector at which the muon entered and exited the detector* The signals 
produced by the detector are electronically amplified and shaped within the 

30 underground detector. The signals are then transmitted to time-to-digitai 
convertors (hereinafter "TDC") units that record the relative times at, and the 
particular wires in which, the signals of each set were produced. 

The detectors employed to ascertain the frequency and trajectory 
of incident muons operate in a manner somewhat similar to a proportional 

35 particle counter. In one particular embodiment of the detectors f the present 
invention, the detector includes a circular, elongate housing that encases a 
plurality of coaxially disposed filter rings that form a plurality of annularly 
shaped gas-containing chambers- The passage of a muon through the chambers 
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results in the liberation of free electrons from the gas molecules. A plurality of 
sensing wires is arranged in a circular array to extend longitudinally through each 
of the chambers* A constant level direct current electrical potential is applied 
between the sensing wires and the corresponding filter rings so that electrons 

5 generated by the passage of a muon through the gas chambers drift at constant 
velocities toward the sensing wires to produce an electrical signal in the wires. 
The relative times at which the signals are produced in the sensing wires are 
indicative of the distances traveled by the electrons and thus reflect the 
distances separating the sensing wires and the path of the muon* These relative 

10 times are monitored by the TDC units for subsequent use by a data processor 
together with the determination of the particular wires in which the signals were 
produced in determining the azimuth angle component of the trajectory of the 
muon. 

To determine the zenith angle and the elevation of the muon 

15 trajectory, a belt composed of electrically conductive pickup wires is wrapped 
around the exterior of the outermost filter ring. The electrical signals produced 
in the sensing wires induce corresponding signals at adjacent locations on the 
pickup wires. The relative times at which the electrical signals reach the 
opposite ends of the pickup wires and knowledge of the precise wires in which the 

20 signals were induced are used to calculate elevations at which the muon entered 
and left the detector. From this information the zenith angle and the elevation 
of the muon may be determined. TDC units are employed to monitor the relative 
times at which the signals reach the ends of the pickup wires- From the zenith 
and azimuth angle components and from knowing the elevation along the 

25 detector at which the muon passed through the detector, the trajectory of the 
muons passing through the detectors may be determined. 

The signals produced on the sensing wires are transmitted to a 
triggering device disposed within the detector, which device makes an initial 
determination of whether the signals possibly originated from a high-energy 

30 muon. 

It is known that the high-energy muons are not affected by 
Coulomb scattering to the same degree as low-energy muons and so high-energy 
muons are deflected relatively little by their collisions with other particles and 
thus typically travel in substantially straight lines. However, low-energy muons 
35 often change direction upon colliding with other particles and, as such, typically 
move along relatively nonlinear paths. Since only the routes which the high- 
energy muons travel along to reach the detector are determinable, only the 
signals produced by high-energy muons are of interest. When the triggering 
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device determines that an electrical signal resulted most likely irom a high- 
energy muon, the TDC units are momentarily stopped and a microprocessor is 
activated to retrieve data from the TDC units. This information along with 
information regarding the location of the detectors below the earth is 

5 transmitted to and stored in a data logger located at the earth's surface. 

As briefly mentioned above, tomographic analysis is used to 
convert data regarding the trajectory and frequency of the muons into the local 
density features of very large underground earth formations. A data processor is 
employed to first determine the generalized or integrated densities of the earth 

10 along incremental angularly separated lines of sight extending through vertically 
'spaced-apart elevations along the detectors and then transform this information 
into localized density values for each spatial point composing the underground 
region under investigation. To reduce the number of computations that must be 
performed 3 the angularly separated lines of sight may in fact be composed of 

15 cone-shaped volumes of sight having their apexes at a particular elevation along 
the detectors and radiating outwardly from such common apex point- Likewise, 
rather than determining a localized density value at every point of the 
underground region, such values can be determined for individual volumetric 
sections of the region. As a result, by the present invention it is possible to 

20 determine the location and size of rather small voids, structural anomalies and 
other density-related features of the earth that may correspond to the locations 
of mineral deposits, promising oil pr natural gas reservoirs or coal horizons. 

Although by the above-described invention it would be possible to 
determine the total density features of the entire volume of the earth through 

25 which the muons pass prior to reaching the detector, this may be somewhat 
impractical since, for a given depth as the zenith angle increases, muons must 
travel greater distances through the earth to reach the detectors. As a result, 
for a particular depth the rate at which the muons strike the detector decreases 
with increased zenith angle. Correspondingly, the time required to gather 

30 sufficient data regarding the density features of the earth surrounding the 
detectors increases with the zenith angle of the muons. So, as a matter of 
practicality, the envelope of the volume of earth investigated by the present 
invention is in the form of a cone having its apex at the location of the 
lowermost detector and its base at the earth's surface. The particular cone angle 

35 that is selected may depend on various factors such as the depth of the 
detectors, the distance to adjacent detection bore locations, and the available 
monitoring time. 
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Brief Description of the Drawings 
The details of one typical embodiment of the present invention will 
be described in conjunction with the accompanying drawings in which: 

FIGURE 1 is a schematic view of a system for investigating the 
5 density features of underground earth formations by monitoring the flux of 
cosmic ray muons penetrating the earth along different trajectories with a muon 
detector shown in isometeric viewj 

FIGURE 2 is a cross-sectional view of the muon detector 
illustrated in FIGURE I; and 
10 FIGURE 3 is a schematic view of a plurality of muon detectors 

constructed according to the present invention and illustrated as disposed within 
a bore hole. 

Best Mode of the Invention 
Referring initially to FIGURE 1, a system 10 constructed according 

15 to the best mode of the present invention currently known to the applicants for 
using cosmic ray muons to ascertain the local density features of large 
underground earth formations is schematically illustrated. System 10 monitors 
the frequency with which the muons pass through a conically shaped section of 
the earth and the trajectory of these muons and then utilizes this information to 

20 develop a three-dimensional map or other representation of the local density 
features of the earth formation under investigation. Additionally referring to 
FIGURE 3, in basic form system 10 includes a plurality of cosmic ray detectors 
12 suspended within a bore hole 14 in interconnected relationship by a support 
line 16. Each time a muon penetrates the earth and passes through a detector 

25 12, a set of electrical signals is produced by the detector that corresponds to the 
presence and trajectory of the muon. The signals are composed of two sets, a 
first set indicative of the azimuth angle of the muon trajectory and a second set 
indicative of both the zenith angle of the muon trajectory and of the elevation 
along the height of the detector at which the muon enters and exits the detector. 

30 The signals are amplified by amplifiers IS and shaped as desirable by signal 
shapers 2k using well known techniques. The amplified and shaped signals are 
transmitted to time-to-digital convertor (hereafter "TDC") units 30 that record 
the relative times at which the individual signals of each set are produced. 
Signals from amplifiers IS and shapers 24 are also transmitted to a triggering 

35 device 36 that determines whether the electrical signals produced by detectors 
12 possibly resulted from a high-energy muon. If this appears to be the case, the 
triggering device momentarily stops the TDC units and activates a microproces- 
sor 38 to retrieve from the TDC units the relative times at which the individual 
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electrical signals were produced by the detectors. This information, together 
with information regarding the depth of the detectors below the ground, is stored 
in a data logger 40 located at the earth's surface and then is transmitted to a 
data processor 42 that utilizes tomographic techniques to convert the informa- 

5 tion stored in the data logger into local density profiles of the underground earth' 
formation. The density profiles may be shown on a visual display unit 44 from 
which geological maps may be made. 

Considering the above, briefly described, method and apparatus of 
the present invention in more detail, as most clearly illustrated in FIGURES 1 

10 through 3, detectors 12 are preferably formed in the shape of an elongate 
cylinder having an outer casing 49 with a tubular body portion 50 and closed ends 
52 to thereby pass down a well or bore hole. Appropriate interconnecting 
devices, not shown, are provided at each end of casing 49 to interconnect the 
detectors with support line and signal bus 16. 

Each detector 12 operates on a principle somewhat similar to the 
operation of a proportional particle counter and, in one preferred embodiment, is 
formed with a cylindrically shaped outer filter ring 60, an intermediate filter 
ring 62 disposed within the outer filter ring, and an inner filter ring 64 disposed 
within the intermediate filter ring. The three filter rings are aligned in coaxial 

20 relationship to each other and extend substantially the full length of casing 49. 
The outer filter ring 60 is smaller in diameter than casing 49 to define an 
annularly shaped outer chamber 66 therebetween, Outer ring 60 and inter- 
mediate ring 62 together define an annularly shaped intermediate chamber 6S 
therebetween, and intermediate filter ring 62 and inner filter ring 64 together 

25 define an annularly shaped, inner chamber 70 therebetween. The chambers 66, 
68 arid 70 are sealed and each is filled with a suitable gas, such as an 
argon/ethane mixture. The specific type of gas may be varied to alter the 
sensitivity of detector 12 in a manner known in the art. 

A plurality of very small diameter sensing wires 72 extends 

30 longitudinally through outer chamber 66. Corresponding sensing wires 74 extent 
longitudinally through intermediate chamber 6g and sensing wires 76exten< 
I 'ongitu -inaHy through intermediate chamber . - 6S and censing wiras 76 exto f* 
longitudinally through inner chamber 70. Except for the number of sensing wire: 
72, 74 and 76 in each chamber 66, 68 and 70, respectively, the arrangement am 

35 construction of the sensing wires are identical; thus, only sensing wires 72 will bi 
described with particularity,' with it being understood that the description of th< 
construction and operation of sensing wires 72 is also applicable to sensing wire 
74 and 76. 
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Sensing wires 72 function similarly to the inner wire electrode of "a 
gas particle counter. They are positioned around the circumference of chamber 
66 and are equally spaced apart from each other to define a circular array. 
Preferably, sensing wires 72 are located radially approximately midway between 
5 casing 49 and outer filter ring 60 and are electrically insulated from casing 49 
and rings 60, 62, and 64. Also, each of the wires 72 is insulated relative to each 
other . 

For detector 12 to operate accurately, it is important that the 
spacing remains constant throughout the entire length of the wires. Thus, if 
10 required, spacers, not shown, may be provided to assist in maintaining the 
spacing between the wires. Ideally the spacers, if utilized, are formed from 
electrically nonconductive material to avoid interfering with the operation of 
detector 12. 

In the operation Of detector 12, an electrical potential is imposed 
15 between wires 72 and casing 49 and outer ring 60 with the wires serving as an 
anode and the casing and outer ring serving as a cathode. Likewise, an electrical 
potential is applied between intermediate sensing wires 74 and outer and 
intermediate filter rings 60 and 62 and between sensing wires 76 and 
intermediate and inner filter rings 62 and 64. The strength of the field in the 
20 vicinity of wires 72, 74 and 76 may be altered by varying the diameter of the 
sensing wires or by varying the applied voltage. As a nonlimiting example, the 
diameter of wires 72, 74 and 76 may be approximately 50 microns. Also, 
preferably the sensing wires are formed from a highly conductive, high-strength 
material, such as stainless steel, but other materials or combinations of 
25 materials may be utilized, such as gold-plated tungsten. 

Preferably, the detector Is operated at a voltage slightly above its 
proportional region in a mode which has been described as a limited streamer 
mode by Atac and Tollestrup in an article entitled " Self -Quenching Streamers" 
Fermi National Accelerator Laboratory, FN-339, 2562.000, (3uly 10, 1981). 
30 When operating in this mode, a large current (approaching 10 milliamperes) of 
relatively short duration can be produced at a specific location in a very short 
time (less than 100 nanoseconds) on the sensing wires. As discussed in this 
article, operation in the limited streamer mode is achieved by proper choice of 
applied voltage, gas mixture and gas pressure. As a typical but nonlimiting 
35 example, the gas within chambers 66, 68 and .70 may be composed of 
approximately fifty percent argon and approximately fifty percent ethane to 
cause the electrons to collect at , a particular location on the sensing wires 
closely adjacent the elevation at which the muon passed by the sensing wires. 
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Also, th gas mixture within the chambers may be pressurized, for instance, to 
approximately two atmospheres to produce electrical pulses on the sensing wires 
that are larger in charge and shorter in duration than if a lower pressure is used. 

When a muon passes through detector 12 along a trajectory, such as 
5 trajectory T illustrated in FIGURE 2, a wake ol free electrons is left behind in 
the gas In chambers 66 and 68- These electrons are attracted to nearby sensing 
wires 72 and 74 by the positive voltage existing in the wires. The electrons drift 
along radial paths toward the adjacent sensing wires, as figuratively illustrated in 
the dot-dash lines in FIGURE 2. As the electrons closely approach the wires, 
10 they gain energy from the electrical field and more vigorously collide with other 
gas atoms and ion pairs causing a multiplication of free electrons. The electrons 
collect at discrete individual locations on the sensing wires 72 and 74 resulting in 
the generation of an electrical pulse or "hit" on each wire, which travels along 
the wire. As discussed more fully below, the existence of the pulse on each 
15 sensing* wire 72 and 74 is independently recorded by a TDC unit 30, as is the time 
at which the pulse occurred. 

The velocity at which the free electrons Initially generated within 
chambers 66 and 68 by a passing muon drift toward sensing wires 72 and 7* is 
dependent upon the strength of the electrical field, the pressure of the gas, and 
20 the type of gas. As a consequence, once these parameters have been chosen, the 
time that it takes for the free electrons to travel through the gas and reach the 
individual sensing wires 72 and 74 is a measure of how close the muon passed by a 
particular sensing wire. This travel time for the electrons can be determined 
from knowledge of the initial time, t Q , at which the muon passed through the 
25 detector and the subsequent times at which electrical pulses were produced on 
the sensing wires. The initial time for the passage of the muon, t n , can be 
calculated from the geometry of the hits that occurred (the particular sensing 
wires in which electrical pulses were produced) and the recorded times of the 
pulses. As long as three or more sensing wires were activated by the passing 
30 muon, the geometry of the hits and the recorded times of the hits can be utilized 
in a "least squares" computation to determine the initial time, t n , for the passage 
of the muon. This type of computation is well known to those skilled in the art 
of the present Invention, however, to the knowledge of the applicants this is the 
first time that this analysis has been carried out in relationship to the particular 
35 geometry defined by the hits occurring on sensing wires 72, 74 and 76 of the 
present invention. 

From the determination f the time interval required for the 
electrons to reach th sensing wires, the distance between the path of the 
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passing muon and the sensing wires can be calculated. With this information and 
with knowledge oi the relative positions of the sensing wires, the precise muon 
trajectory projected on a transverse plane or azimuth angle, can be determined. 

The zenith or dip angle component of the muon trajectory and the 
5 elevation at which the muon enters and exits the detector (hereafter "elevation") 
are also sensed by detector 12. With the known zenith and azimuth angles and 
elevation of the muon, the trajectory of the muon can be determined. The 
components of detector 12 that are used to determine the zenith angle and 
elevation of the muons include a belt 100 of cathode pickup wires 102, which is 
10 wound in spiral fashion around outer filter ring 60 to form convolutions or 
windings which cover the outer surface of the filter ring* The ends of the belt 
may be beveled to substantially completely face the upper and lower end portions 
of ring 60. The pickup. wires 102 comprising belt 100 are equally laterally spaced 
apart from each other, disposed parallel to each other and extend the entire 
15 length of the belt having opposite ends at the upper and lower ends of the 
detector. The wires are electrically insulated from each other and from filter 
ring 60, As with sensing wires 72, 7fc and 76, pickup wires 102 are preferably 
constructed from durable, highly conductive material, such as copper. The upper 
and lower ends of each wire 102 are connected to a TDC unit 30 through the 
20 inter mediacy of an amplifier 18 and a shaper 24. 

In operation, when a muon passes through detector 12 and causes 
electrons liberated from the gas in chamber 66 to produce an electrical pulse in 
sensing wires 72, corresponding pulses are induced in pickup wires 102 at 
locations adjacent the elevation in which the electrical pulses are produced in 
25 the sensing wires. The electrical signal induced in pickup wires 102 travels 
upwardly and downwardly along the length of the pickup wires. By using the TDC 
units 30 to monitor the relative times at which the electrical pulse reaches the 
opposite ends of each wire 102, it is possible to determine the particular 
convolution or wrap at which the pulse was induced in the pickup wires. Also, 
30 from the wrap angle of the belt 100, It is possible to compute the location along 
the length of the convolution at which the pickup wire is located most closely to 
the particular sensing wire 72 from which the electrical pulse was induced in the 
pickup wire. With this information, it is possible to calculate the precise 
location along the length of pickup wire 102 at which the electrical pulse was 
35 induced and thus determine the elevations along the height of detector 12 at 
which the muon entered and exited the detector. This in turn enables the zenith 
angle and the elevation to b ascertained. 

As an alternative to sensing the difference in time that the 
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electrical signals reach the ends of pickup wires 102 lor determining the zenith 
angle and elevation of the muon, instead the size of the electrical charges 
reaching the ends of the signal wires can be measured. The relative size of the 
charges reaching the opposite ends of the wires is indicative of where along the 
5 length of the wire the hit occurred. From this information, the particular 
convolution or wrap at which the pulses were included in the wire can be 
determined. 

The number of pickup wires 102 comprising belt 100 may be varied 
to optimize the operation of detector 12. As a muon passes through chamber 66, 

10 electrical signals may likely be produced on a cluster of sensing wires 72, and the 
electrical pulses produced on each sensing wire 72 may induce a corresponding 
electrical pulse or hit on pickup wires 102. In wider belts composed of a larger 
number of wires, there is less likelihood of more than one hit occurring in a 
particular pickup wire and thus less likelihood of miscalculating from what 

15 particular sensing wire the electrical signal in the pickup wire was induced. 
However, because with each convolution of the belt each pickup wire 102 extends 
a longer distance along the length of the detector in wider belts than in narrower 
belts, the particular location along the length of the detector at which the signal 
was induced in a pickup wire is less accurately determinable. Also, a larger 

20 number of amplifiers 1&, shapers 2k and TDC units 30 are required. In the 
opposite situation, if a narrow belt of wires 102 is used or if belt 100 is composed 
of a single pickup wire, more than one hit will occur in a particular pickup wire, 
thereby creating ambiguities as to what particular sensing wire 72 induced the 
hit in the pickup wire. 

25 Rather than being wound around outer filter ring 60 in the manner 

illustrated in FIGURES 1 and 2, belt 100 instead may be disposed in other 
locations within detector 12 such as along the inside surface of casing 1^ or 
around the outside surface of intermediate filter ring 62. The elongate, 
cylindrical shape of the above detector 12 is particularly well suited to be placed 

30 down vertical or diagonal wells or bore holes. It is to be appreciated that if a 
detector is located elsewhere, such as in a mine drift or shaft or adjacent a cliff, 
the detector can be constructed in other shapes, such as rectangular, without 
departing from the spirit or essential characteristics of the present invention. 

As discussed above, because low-energy muons are significantly 

35 affected by Coulomb scattering and thus change directions upon colliding with 
other particles, it is n t possible to trace the path through which they passed 
through the earth before reaching the detector. To help reduce the numbers of 
low-energy muons which pass through chambers 66, 68 and 70, to thereby help 
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minimize the frequency at which nonusable signals are generated in sensing wires 
72, 74 and 76 and in pickup wires 102, preferably filter rings 60, 62 and 64 are 
formed from high-density material such as lead. 

System 10 also includes amplifiers IS and shapers 24 located within 

$ casing 49 for amplifying the signals received from sensing wires 72, 74, and 76 
and pickup wires 102 and altering the signals into a shape more appropriate for 
triggering device 36. Such amplifiers and shapers are well known in the art and 
thus will not be further described. 

The amplified and shaped signals are transmitted to TDC units 30 

10 that record the time at which pulses are produced in sensing wires 72, 74 and 76 
and the times at which induced pulses reach the opposite ends of each pickup 
wire 102- The TDC units continually operate unless momentarily stopped to read 
out the information stored within the units, as discussed more fully below. TDC 
units that are capable of recording events within an accuracy of a few 

1 5 nanoseconds are now available. An example of one such TDC unit is 
manufactured by Radio Corporation of America and designated as Model 
No. TCS282. The TDC units are located within casing 49 and are interconnected 
with amplifiers 28 and shapers 24 by transmisison lines. 

The amplified and shaped signals are also transmitted to a 

20 triggering device 36 through the use of a triggering bus 110- Triggering device 
36 makes a determination of whether the signals might have originated from a 
high-energy muon. If the signals could have resulted from a high-energy muon, 
the triggering device activates a microprocessor 38 that momentarily stops the 
operation of the TDC units so that the information stored within the units may 

25 be read out by the microprocessor. After this information is received from the 
TDC units, the microprocessor permits the TDC units to return to their normal 
operation to again record the times at which electrical signals are produced in 
sensing wires 72, 74 and 76 and in pickup wires 102. As a result, the TDC units 
are permitted to operate proportionally more of the time in their data intake or 

3q record mode than in their readout mode. This is important since, during the 
times that the microprocessor is retrieving information stored in the TDC units, 
the TDC units cannot record hits occurring on the sensing or pickup wires. 
Moreover, by using the triggering device to attempt to discriminate between 
signals caused by high-energy muons and other sources, the proportion of 

3 5 nonusable data reaching data logger 40 and data processor 42 is reduced. 

To enable it to function in the manner described above, triggering 
device 36 may be designed in various ways. For instance, the triggering device 
may be designed to activate microprocessor 38 when signals are produced on 
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more than one sensing or pickup wire. Alternatively, the triggering device may 
be constructed to activate the microprocessor when electrical signals ar 
produced on the wires in a particular pattern that could be equivalent to making 
an initial determination of the trajectory of the muon. 

5 To reduce the complexity of the triggering device, it may be 

designed to distinguish signals produced from high- and low-energy muons from 
signals caused by background radiation and noise since this is an easier 
determination to make than trying to distinguish between the signals caused by 
high- and low-energy muons. In this situation, the high- and low-energy muons 

10 are differentiated from each other when their trajectories are calculated by data 
" processor 42. 

The minimum number of muons that must pass through detectors 12 
from a given direction for system 10 to be able to determine the mean density of 
the underground earth formation along the individual angularly separated lines of 

15 sight depends on the desired level of accuracy of the density determination. 
Factors that affect the rate at which muons may strike a* detector 12 from a 
particular angle include the density and composition of the surrounding 
underground earth formation and the distance that the muons must travel 
through the earth to reach detectors 12- This distance depends on, inter alia, the 

20 depth of the detectors below ground level, the topography at the surface of the 
earth and the zenith angle of the muon. Although muons strike the earth*s 
surface from all angles and not just from the radial direction, as the zenith angle 
of the muons increases, the distance through the earth which they must travel to 
reach the detector 12 also increases. As a result, for a particular detector depth 

25 the number of muons reaching detector 12 normally diminishes with an increasing 
zenith angle. Thus, as a matter of practical consequence, preferably the volume 
of the earth that is investigated by system 10 is shaped in the form of a cone 
having its apex located at the bottom of the lowermost detector 12 and its base 
located at the surface of the earth. The cone-angle of the cone may be in the 

30 order of approximately 60 degrees. However, the particular cone-angle chosen 
will depend on such factors as the depth of the detectors and the time available 
for monitoring the muons. It is to be appreciated that the present invention 
enables volumes of earth to be investigated that are several orders of magnitude 
larger than can be investigated by currently known methods, including the use of 

35 gamma ray detectors. 

As briefly discussed above, triggering device 36 signals the 
microprocessor when a possible high-energy muon has passed through detector 
12. When notified, microprocessor 38 momentarily stops the TDC units 30 so 
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that the latest data monitored by the TDC units can be extracted and 
transmitted to a data logger 40 located aboveground. By retrieving only the 
information received by the TDC units in the latest time interval of a 
predetermined length, extraneous information from earlier arriving low-energy 
5 muons, background radiation or other noise sources is not retrieved by the 
microprocessor. As with triggering device 36, microprocessor 38 is also located 
within detector casing 49, The microprocessor is interconnected with a data 
logger 40 by an appropriate transmission line in a manner well known in the art. 

Since muons do not strike the earth with sufficient frequency for 
10 density determinations to be made in real time, a data logger 40 is used to store 
information received from microprocessor 38 for subsequent analysis. Data 
loggers, such as logger 40, are in common use and thus need not be further 
described. They typically store information on different types of hardware, such 
as on a magnetic tape or disc. 
15 System 10 also includes a data processor 42 that uses the data 

stored in logger 40 to determine the trajectory and flux of the muons passing 
through detectors 12. As discussed above, the muon flux reaching a detector 
along a particular trajectory depends upon the density of and the distance 
through the earth that the muon travels through to reach the detector. Thus, 
20 from knowledge of the trajectory of the muons passing through the detector, the 
flux at each of these trajectories and the distance through the earth that the 
muons travelled, the data processor can determine integrated density values 
along individual lines of sight passing through a specific elevation or elevational 
range of the detector. Each line of sight has an azimuth angle component and a 
25 zenith angle component and divides the underground region into discrete paths 
that pass through the particular elevation or elevational range of the detector. 
Alternatively, the lines of sight may in fact be in the form of individual cone- 
shaped volumes having their apexes at the detector and radiating outwardly 
through the earth formation from the detector. The integrated density values 
30 lor each of these lines of sight is the result of the totality of the densities of the 
different portions of the earth formation through which the lines of sight extend. 
These integrated density determinations along the lines of. sight that the 
underground formation is divided into may be made at spaced-apart elevations 
along the length of each detector or along discrete elevational ranges of each of 
35 the det ctors. By determining integrated density values at different points along 
the detector(s), muons passing through a particular point or location in the 
underground region may be monitored at different angels in a manner similar to 
which a tissue structure in a human body is viewed from many different X-ray 
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angles during a computerized axial tomography (hereafter "CAT") scanning 
examination. In the present situation, the data processor can be programmed to 
perform a tomographic analysis to thereby combine the totality of the integrated 
density values determined for the lines of sight for all of the elevations or 

5 elevational ranges to thereby generate the localized density values for each 
specific spatial point in the underground region. To reduce the numerical 
computations which must be made and/or to reduce the amount of data which 
must be accumulated, rather than determining the local density value for each 
specific point a density value can be determined for discrete volumetric sections 

10 of the underground region. As in CAT scanning, these individual density values 
' can be unfolded from the integrated density values that represent a summation 
of/or the total density of the earth along each of the lines of sight. 

By using the method and apparatus of the present invention to 
determine localized density values for each individual point or volumetric 

15 sections of the region being investigated, different density features of the region 
may be distinguished from each other, even those that are closely located 
together and/or are approximately of the same density. Also, the positions of 
the density features relative to the detector may be determined. In general, the 
use of tomographic analysis makes it possible to accurately and efficiently 

20 detect localized density related features of very large volumes of earth that 
heretofore has not been possible. Due to the precision and exactness with which 
density related features may be investigated with the present invention, it may 
serve as a valuable tool for the exploration of petroleum reserves, mineral 
deposits, water or other subjects of interest that exhibit characteristic density 

25 features. 

The density features of the underground as generated by data 
processor 42 may be displayed in various quasi three-dimensional forms on visual 
display 44, such as in a transaxial image format or a generalized planar image 
format. In addition, the results of the tomographic analysis can be transmitted 

30 to a printing mechanism to generate geological maps of the underground region 
being investigated. 

As discussed above, the distance that muons travel through the 
earth is one of the factors utilized in determining the integrated density values 
for the angular sights radiating, from each detector. Because the earth surface is 

35 not perfectly smooth, the distances travelled by the muons to reach a detector 
are depefident upon the topography of the earth in the area surrounding the bore 
hole. Thus, it is important that this topography is studied and the results of the 
study incorporated into the mean density value calculations. 



-19- 



0081314 



It is to be understood that data processor 42 need not be located at 
the site of bore hole 14. Also, the tomographic analysis does not have to occur 
concurrently with the arrival of information at data logger 40, rather it can take 
place at a later time, lor instance when sufficient data has been gathered or 

5 when computer facilities are available* 

As will be apparent to those skilled in the art to which the 
invention is addressed, the present invention may be embodied in forms and in 
methods other than those specifically disclosed above without departing from the 
spirit or essential characteristics of the invention. The particular steps lor 

10 ascertaining density features of underground earth formations and the particular 
embodiments of detection system 10, as described above, are therefore to be 
considered in all respects as illustrative and not restrictive, with the scope of the 
present invention being set forth in the appended claims rather than being 
limited to the foregoing description. 



-20- 



0081314 



CLAIMS 

L A method for ascertaining density features of underground 
earth regions, comprising: 

placing a muon detector at an underground location; 
5 sensing the frequency and trajectory of incident muons striking the 

detector; and 

determining the density features of the underground regions based 
on the frequencies and trajectories of the muons striking the detector. 

2. The method of Claim 1, wherein the step of determining the 
10 density features of the underground region includes determining integrated 

density values for incremental, angular paths along which muons pass through 
specific elevations along said detector. 

3. The method of Claim 1, wherein the step of determining the 
density features of the underground region further includes combining all of the 

1 5 integrated density values by tomographic analysis to generate density values for 
individual spatial locations of the underground region being investigated. 

% m The method of Claim 1, wherein the step of sensing the 
frequencies and trajectories of incident muons striking the detector includes 
producing output signals corresponding to the presence of each muon and its 
20 trajectory; and the step of establishing the density features of underground 
regions includes processing the output signals. 

5. The method of Claim I, wherein the step of sensing the 
frequencies and trajectories of the muons includes producing a first set of 
electrical signals corresponding to a first angular coordinate of a muon 
25 trajectory and then utilizing the first set of electrical signals to induce a second 
set of electrical signals corresponding to a Vertical coordinate and to a second 
angular coordinate of the muon trajectory. 

6- A detector module for detecting the flux and trajectories of 
muons, comprising: * 
30 a chamber containing ionization gas; 

sensing wires disposed within the chamber and arranged so that * 
when a muon passes through the chamber electrical signals are produced in the 
wire which are indicative of a first vector component of the muon; and 

pickup wires disposed angularly to the sensing wires, in which the 
35 pickup wires' electrical signals are induced by the signals in the sensing wires 
that are indicative of a second vector component f the muon passing through 
the chamber. 
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7. The detector module according to CJaim 6, including a 
plurality of gas chambers disposed adjacent each other, with a plurality of 
sensing wires disposed within each of the gas chambers. 

X. The detector module according to Claim 6, wherein said gas 
5 chamber is annular in shape; and said sensing wires are arranged in a circular 
array to extend longitudinally through said chamber. 

9. The detector module according to Claim 6 or 8, wherein said 
pickup wires are arranged in a belt composed of a plurality of spaced-apart 
pickup wires, said belt at least partially covering said gas chamber. 
10 10. A system for ascertaining the local density features of 

underground earth regions, comprising; 

muon detector means lor producing a plurality of sets of electrical 
signals when a muon passes through said detector means, said signal sets being 
indicative of the trajectory of the muon; 
15 recorder means for recording said signals; 

a microprocessor; 

a triggering device for determining whether the signals produced by 
said detector means possibly originated from a muon or from some extraneous 
source, said triggering device activating said microprocessor to retrieve said 
20 signals from said recorder if said signals possibly originated from a muon; 

data processing means for converting the electrical signals re- 
trieved from said recorder means by said microprocessor into the local density 
features of the underground earth region- 
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